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Abstract 

Effect of long range charge density oscillations on capacity of tunnel 
structure is studied within random phase approximation (RPA). Us- 
ing this approximation we obtain expressions for the short and the 
long distance behavior of the self-consistent screened potential. We 
demonstrate that long range density oscillation, commonly referred to 
as Friedel oscillations, leads to decreasing of total electrostatic capacity. 
Particular emphasis has been placed on influence of an external mag- 
netic field is applied perpendicular to the barrier plane on the capacity 
of a structure. It is shown that increasing of magnetic field implies an 
increase of quantum correction to capacity due to Friedel oscillations. 

1 Introduction 

Response of tunneling structures on an external dc bias plays an important role in 
a number of tunneling phenomena, such as zero-bias anomaly and Coulomb block- 
ade. Besides only the fundamental physical interest, it can have a pronounced 
effect on performances of electronic devices, based on semiconductor tunnel junc- 
tions. 

It is well known that the presence of a potential barrier leads to the long-range 
oscillation of electron density, which are commonly referred to as the Friedel oscil- 
lations. In the contrast to the three-dimensional case where the electron density 

oscillations around an impurity decays as ^'^^^'^^^^) (where r is the distance from 

the impurity and ky is Fermi wave vector) in the case of metal surface this oscil- 
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lations decays only as here z is the distance from the surface [1]. In [2] 

was shown that in one ID tunnehng structure, charge density oscillations decay 
only as '"^^(^^'f^) g^j^j ^j-^jg induces a singularity in differential conductance. This 
singularity originates from the electron scattering on the Friedel oscillations. 
It is well-known that the motions of three-dimensions electrons in an external mag- 
netic field are quasi one-dimensional and it is possible to expect that in a strong 
fields the long-range behaviour of electron density will be similar to those of one- 
dimensional systems. 

In this paper we investigate effect of Friedel oscillations on capacity of barrier 
structure, it will be shown, that long-range electron density oscillations lead to 
decreasing of electrostatic capacity of the structure and this effect plays more im- 
portant role with increasing external magnetic field, which is perpendicular to the 
barrier plane. 

Formally, the goal is to find a self-consistent distribution of a total electrostatic po- 
tential $(r), since the knowledge of this potential is sufficient to calculate capacity 
of a structure. Self-consistent potential is calculated within RPA approximation, 
using this approach we obtained the short- and the long-distance behaviour of 
electrostatic potential. 

2 General formulation 

Here we briefly summarize the standard self-consistent analysis, which allows to 
write an expression for potential $(r) . In linear response theory the induced 
charge density is given by 

/oo 
dr'n(r,r')$(r') (1) 
-oo 

here the polarization kernel n(r, r') is given by 

n(r, r') = E Tr^^*^^9^^k(r')*k'(r')^k.(r)^k(r). (2) 
k,k' -^k' — rjk + tv 

Hence fk is the Fermi distribution function (below we will assume that the tem- 
perature iz zero), is wave functions. 

The potential $(r) is also related to the charge density by Poisson's equation which 
can be written as 

^r) = -ATI J Gir,r')gir'), 
here G{r, r") is a Green function of Laplas equation. 

Combining (1) with (2) gives us an integral equation for electrostatic potential 
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here 



$(r) = J dr'F(r,r')$(r'), 
F(r,r') = J dr"G'(r,r")n(r",r'). 



(3) 



as it was mentioned in introduction, the goal of the present work is to calculate the 
the self-consistent response of tunnel structure on external bias, because formally 
the knowledge of potential is sufficient to calculate capacity of a structure. 

The model wc use for a tunnel structure consists two semiconductor half-spaces 
z < —d/2 and z > d/2, separated by a potential barrier of width d. A uniform 
compensating positive background extending through semiconductors, so electro- 
static potential $(r) doesn't depend on lateral (x-y) coordinates. Below we will 
assume that: 

(i) transparency of the barrier is zero, so electrostatic potential is a linear func- 
tion inside a potential barrier. 

(ii) ^{z) is an odd function in z direction. 

This assumptions allow to transform the task to the problem of an electrical field 
E penetration at the surface of semiconductor [3], [4]. However in the difference 
from the case of surface of semiconductor, where E had a fixed value, in our case 
the electrical field depends on external dc bias V as well as induced electrostatic 
potential. The total electrical field E in the barrier structure can be written as a 
sum of two terms E{z) = E^xtiz) + -Eind(-2), where E^xtiz) is discontinous function 
which is equal to a constant E^xt inside the barrier and zero in semiconductor lead- 
ers, £'ind(-2) is the field induced by electron systems and it is discontinous function 
too. The fact, that a total electrical field should be continues function along the 
structure, bring us to the boundary condition for E±-a_dL at {z = d/2 + 0) 



It was mentioned that self-consistent potential is calculated within RPA approxi- 
mation, which is valid if rg <^ 1. Where 



here is Bohr radius In this study we ignore considerations of exchange and 
correlation interactions between electrons. 

3 Self-consistent potential of tunnel structure 

In this section, we present calculations of the self-consistent potential for a tun- 
neling structure. 




(4) 
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The case of zero magnetic field 

In this case electrons are described by a normalized eigenfunction 

where ky and ry are two-dimensional vectors laying in the x-y plane parallel to 
the barrier layers and z-dependent part of wave function is ipk^ — . — sin k^z. On 
substituting this wave functions in Eq.(3) we obtain 

/•oo 

/ dz'[A{z - z') + A{z + z' - d)]^z')dz'- 
^z) = ■'i^ z> dl2 (5) 

/ K{z,z')^{z')dz' 

Jd/2 

^"{z) =0 1^1 < d/2 (6) 

here 

dk\\ r°° dkdk' f e'^'"'''"'^ + c.c.\ fk',k^^+g - fk,k^^ 



/oo 
-oo 



J-oo ZTT r JO 



oo(27r)2 7o (2^)2 \(k-k'y + q^jEk',k,^+,-Ek,k^^+iO 
°° dk\i ["^ dkdk' fk',k\^+q - fk,k\^ 



/ ^i{k+k')z^-i{k-k')z' _|_ ^i{k+k')z ^i{k-k')z' _|_ ^i{k-k')z ^i{k+k')z' _|_ ^i{k-k')z^-i{k+k')z' ^ ^ ' 

[ {k + k'f + g2 + {k - k'f + g2 ^ 

here q = (g^, qy, 0) is 2D wave vector in the plane of barrier. As far as $ doesn't 
depend on lateral coordinates, we will consider the case q = 0. 
Making Fourier transforming Eq.(5) with respect to the variable and some algebra 
we obtain 

*(C) --TT / dq,^--K{q,,riMr)) ( > 0, (7) 

Zn Jo J-oo Qz^kQz) 

here ( = z — d/2, rj = z' — d/2, e{qz, 0) is the Lindhard dielectric function 

/oo 
dze-"i^'A,{z), 
-oo 

/oo _ 
dr^e-'^^^K,{C,r]). 
-oo 

Integral equation (7) can be considered as a generalization of well-known Shafra- 
nov's equation [5] to the case of non-uniform electron gas. 
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After integrating over and k, we obtain the following expression for the integral 
kernel 



'''^ \Qz\ \ , 1^ ■ ni , sin 2A;f?7 sin 5^77 cos 2k^r] sin q^r] 

[ +2^""'^^^^°^^^^ + ^^^f M ' 

(8) 

here /cxf is Thomse-Fermi wave vector. In the long wavelength limit [q^ "C k-p) 
RPA dielectric function and integral kernel Kcan be simplified as 

<^z) = 1 + ^ 

r^/ N 2 sin(2A;Fr?) 
K{qz,V) = -^k^^F—^i cosg^r/. 

In this limit the integral over in (7) is done by closing the integration contour 
in the upper half plane and picking up the pole at q^ — ikj^ and Eq.(7) can be 
transformed to 



|2 



with boundary conditions 

Web Wcf) 

= (10) 



C=o 



E * 



C=oo 



In deriving of Eq.(9) we used the following expression 

Here we emphasize that under our assumption (i) electron density must be zero at 
interface semiconductor-potential barrier C = and as it easy to see that solution of 
Eq.(9) satisfies this condition. It describes the behaviour of the screened Coulomb 
potential if C ~ ^tf and it exponentially decays at large distances ( ^ Ate- 
Here we would like to repeat that Eq.(9) was obtained with account of pole in e, 
however apart from the pole, the dielectric function has branch point at q^ = ^k^ 
and this singularity generates long-range oscillations in charge density (Priedel 
oscillations). At large distances C ^ Ap this oscillations plays the dominant role in 
distribution of induced charge. Asymptotic solution of Eq.(7) can be presented as 

^ ^sin2^ 
A;f (2A;fC)' 
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Thus behaviour of electrostatic potential across the structure is 



-fcxF / dx\ 1 

$(0)e -^0 V ^hx 



^Q)k^^j^^^0{k^/h) C>Cc 



(11) 



According to [4] parameter C,c can be estimated as 



Zc = ATF(ln^ + a), 



TF 

where a > 1. 

The case of finite magnetic fields 

In this section we consider the response of barrier structure on external bias in 
the presence of magnetic fields perpendicular to the planes of the barrier. Under 
the Landau gauge with vector potential A — {—Hy, 0, 0) the wave functions and 
corresponding energy levels can be specified by the set of quantum numbers {n, kz) 
as ^ 

*a(r) = ^e^'^^iJkMMy - yo), 



En{kz) = h(jJc{n + 1/2) + 
eH ■ 



2m 



Here n is a number of Landau level , oUc — is the cyclotron frequency and 0n(y) 
is the normalized harmonic-oscillator wave function with Landau state index. 
As was mentioned above, applied magnetic field reduces the effective dimension- 
ality of charge from 3D to ID, in other words electron gas in semiconductor can 
be considered as a set of ID gases and every one dimensional gas is specified by 
partial Fermi wave vector k'^. In the case of sufficiently weak fields k^l^ >> 1, the 
partial Fermi wave vector can be defined as 



:i- 



(n+1/2)). 



Using together with the expression for polarization operator and we get an equation 
for self-consistent potential in external magnetic field. 

1 roo POO pilzQ 

27r Jo J-oo qien{qz) 
where eniqz) is dielectric function in the presence of magnetic field. 



i + E 



aBlhz 



In 



2k^^ + qz 
2A;f" - qz 



(13) 
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In the long-wave limit <^ the expression for en can be transformed to 

eH(g.,0) = l-54^, (14) 

here k^"^ = , k1 = J2nks'^ and summation is performed over all Landau 

levels where k^ <C k^. It is evident that this condition isn't valid for all levels, 
however using numerical modeling of (13) it is possible to demonstrate that in 
a case of weak-magntetic fields expression (14) is a qiute good approximation of 
(13), even though the summation is performed over all occupied Landau levels. 
It is easily seen that in the zero-magnetic-field limit kg transforms to well-known 
expression for Thomas-Fermi vector. 

In a similar to the zero-field case, in the long- wave limit the integral equation 
(12) can be transformed to 

- ('^s' - E cos 2kiX)^ = 0. (15) 



E — 



(16) 



C=o 

This equation describes behaviour of electrostatic potential at C ~ Ag, here is 
Thomas-Fermi length in magnetic field. It is possible to show (using Poisson sum- 
mation formula) that in the limit of zero magnetic field Eq.(15) can be transformed 
to Eq.(9). 

Similarly to the case H = 0, long-range behaviour of $ is conditioned by branch 
points in the dielectric function and at C ^ Atf charge density oscillations have 
the following asymptotic 

g{c) = mEk:kf-^l^. (17) 

Every term in this sum decays as 1/C and it is typical for one-dimensional sys- 
tems. The same behaviour was obtained in [6], however in that article authors 
used exponential parametrization of self-consistent potential, where unknown pa- 
rameter determined while functional minimization process. In the case l^kp ^ 1, 
summation in (17) can be fulfilled using Poisson formula 

The first term in this expression corresponds to density oscillations without mag- 
netic field and it is in an agreement with Eq.(ll), the second term depends on 
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magnetic field. 

In the magnetic quantum limit, when only the lowest Landau level is partially 
filled, we have 

.(c) = mk.k.'^^ (19) 

here we need to mention that in quantum limit, can't be calculated using ex- 
pression presented above. 



4 Electrostatic capacity of the barrier structure 

It is well-known that the differential capacitance per unit area is C = dQ/dV 
where _ 

poo 

Q ^ dzg{z) 

Jd/2 

is a total charge over semiconductor lead. Substituting distribution of self-consistent 
electron density g{z) into the definition for electrostatic capacity we will obtain an 
expression for C. In the case of zero magnetic fields it expressed as 

^"4^(d + 2ATF + 2A,)' 



Aq = rsln ^ ^= Atf- 



where 



Accounting of Priedel oscillations in electron density produce a shift of charge cen- 
ter mass towards and away from the surface semiconductor-dielectric, what leads to 
decreasing capacity of the structure. We wish to separate the classical term in the 
total capacity, which depends on the barrier width and classical screening length 
from quantum one which depends on long-range oscillations of charge density 

1111 



z' ' 
^ *-^q 
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here Cg is the geometrical capacity of the structure, Cg = -— - — is contribution to 

capacity due to Thomas-Fermi screening in semiconductor regions. The expression 
for Cg coincides with result obtained within semiclassical approach [7], which was 
based on solution of Boltzmann equation for electron distribution function. Cq = 

is quantum correction to tunnel capacity due to long-range behaviour of 

87rA„ 
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charge density. As we mentioned in the introduction, in this study we assumed 
that Ts <^ 1, so quantum contribution is smaller then the classical one 



^ = ^ « 1- (21) 



In the case of finite magnetic fields Friedel oscillations decay slower than in the 
case H — and this fact brings to increasing to penetration length of electric field 
in semiconductor lead. For this case the expression for capacity has the form 

^'^i ^'2) 

n kg 

Under the condition kpln ^ 1 we have 

Ch = 1^ ^TTfc (23) 

d + 2Atf + Aq + 2ATF7f^ E sinnkihly^f 

Thomas-Fermi screening length also depends on H, but in a weak magnetic fields 
it is approximately equal to Atf. 



Cs 

As it is easy to see that — — — - is still small however in difference from the previous 

case this ratio ~ y/rs- Quantum contribution in the total capacity will play more 
significant role with increasing of magnetic field. 



5 Conclusions 

Within RPA approximation expressions for self-consistent electrostatic potential 
are obtained. It is shown, that long-range electron density oscillations lead to 
decreasing of structure capacity and this effect plays more significant role with 
increasing external magnetic field. 
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